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1. INTRODUCTION 



Consider the matrix Schrddinger equation on the half line 

-ip" + V(x) ip = k 2 i(j, iGR + , (1.1) 

where R + := (0, +00), the prime denotes the derivative with respect to the spatial co- 
ordinate x, and the potential V is a n x n selfadjoint matrix-valued function integrable 
in x. The integrability V G L 1 (R + ) means that each entry of the matrix V is Lebesgue 
measurable on R + and 

/■OO 

/ dx\\V{x)\\ < +00, (1.2) 
Jo 

where ||V(a;)|| denotes a matrix norm. Since all matrix norms are equivalent, without loss 
of generality, we can use the matrix norm defined as 

n 

\\V(x)\\ :=max^|y ;s (x)|, Z = l,...,n, 

s=l 

where Vi s (x) denotes the (/, s)-entry of the matrix V(x). Clearly, a matrix- valued function 
is integrable in x if and only if each entry of that matrix belongs to L 1 (R + ). 

Note that V is not assumed to be real valued but is assumed to be selfadjoint, i.e. 

V(x) = V(x)\ x G R + , (1.3) 

where the dagger denotes the matrix adjoint (complex conjugate and matrix transpose). 
The wavefunction tp(k, x) appearing in (1.1) will be either an n x n matrix-valued function 
or it will be a column vector with n components. 

When the analysis of (1.1) is considered at or near k = 0, in addition to the integra- 
bility we require that the potential V has a first moment, i.e. 

POO 

/ dxx\\V(x)\\ < +00. (1.4) 
Jo 
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All the results in our paper are valid under the assumption that the potential V is selfad- 
joint and belongs to L}(R + ), i.e. 

/•OO 

/ dx(l + x)\\V(x)\\ < +00. (1.5) 
Jo 

We are interested in studying (1.1) with a selfadjoint potential V in L\(R + ) under 
the general selfadjoint boundary condition at x = 0. As indicated in [2], without loss of 
generality, it is convenient to state the general selfadjoint boundary condition at x = for 
(1.1) in terms of constant n x n matrices A and B such that 

-BV(0) + ^V(0)=0, (1-6) 

-B t A + A t B = 0, (1.7) 
A ] A + B^B > 0, (1.8) 

i.e. A^ B is selfadjoint and the selfadjoint matrix (A^A+B^B) is positive. There are various 
equivalent formulations [2,10-14] of the general selfadjoint boundary condition at x = 0. 
Let us also mention that it is possible to use some transformations on A and B without 
affecting (1.7) and (1.8). In Section 4 we will elaborate on two such transformations, 
namely (A, B) (AT, BT), a right multiplication by an invertible matrix T, and a unitary 
transformation (A, B) h-> (M^ AM, M^BM), where M is a unitary matrix. A combination 
of such transformations as given in (4.10) will turn out to be useful. 

We have two primary goals in this paper for the Schrodinger equation (1.1) with the 
selfadjoint boundary condition (1.6)-(1.8) when the potential V satisfies (1.3) and (1.5). 
Our first primary goal is, even when V satisfies the weaker condition (1.2) instead of 
(1.5), to establish the large-/c asymptotics of various quantities related to (1.1) such as 
some relevant scattering solutions, the Jost matrix, the inverse of the Jost matrix, and the 
scattering matrix. Our second primary goal is, under the additional assumption (1.4), to 
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derive so-called Levinson's theorem, namely to obtain the relationship between the number 
of bound states and the change in the phase of the determinant of the scattering matrix. 

A bound state corresponds to a square-integrable column- vector solution to (1.1) satis- 
fying the boundary condition (1.6)-(1.8). The selfadjointness (1.3) of V and the selfadjoint 
boundary condition (1.6)-(1.8) assure that the corresponding Schrodinger operator is self- 
adjoint on L 2 (R + ), and hence its eigenvalues must be real. When k 2 > 0, it turns out 
that there are no square-integrable column-vector solutions to (1.1). As a result, a bound 
state, if it exists, occurs only when k 2 is negative, or equivalently when k appearing in 
(1.1) is on the positive imaginary axis in the complex plane C. Thanks to the restriction 
(1.5), the number of such /c-values turns out to be finite. We will see that, at each /c-value 
corresponding to a bound state, the number of linearly independent square-integrable 
column-vector solutions (i.e. the multiplicity of the corresponding bound state) cannot 
exceed n. The number of bound states is defined as the number of bound states including 
the multiplicities. 

The large-/c analysis and Levinson's theorem for (1.1) with the selfadjoint boundary 
condition (1.6)-(1.8) are relevant in the study of the corresponding direct and inverse 
scattering problems. The direct scattering problem for (1.1) is to determine the scattering 
matrix and the bound-state information when the matrix potential V and the selfadjoint 
boundary condition are known. On the other hand, the inverse scattering problem is to 
recover the potential and the boundary condition from an appropriate set of scattering 
data. 

Our paper complements the study by Agranovich and Marchenko [1] , where the large- 
k asymptotics and Levinson's theorem are provided only under the Dirichlet boundary 
condition. Our study also complements the study by Harmer [10-12], where the general 
selfadjoint boundary condition is used but the large-/e asymptotics of the scattering matrix 
is obtained by only providing the leading term with the remaining terms as specified as o(l) 
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as k becomes large. In our paper, we not only provide the leading term but we also specify 
the next order term and establish the large-/c asymptotics up to 0(1/ k 2 ), which is crucial 
in establishing the Fourier transforms of various quantities relevant to the corresponding 
inverse scattering problem. 

Our current paper also complements our own recent study [2], where the rigorous 
small-/c analysis for (1.1) is provided with the general selfadjoint boundary condition (1.6)- 
(1.8). The small-/c analysis in [2] is crucial in the derivation of Levinson's theorem in our 
current paper. In fact, the only reason for needing (1.5) rather than merely (1.2) in 
our current paper is because of the fact that the small-/c asymptotics are also needed to 
establish Levinson's theorem, and those asymptotics require (1.5). 

The half-line matrix Schrodinger equation (1.1) has applications in quantum mechan- 
ical scattering involving particles of internal structures such as spins, in scattering on 
graphs [4-9,15-18], and in quantum wires [13,14]. The consideration of the general selfad- 
joint boundary condition at x = given in (1.6)-(1.8) rather than the Dirichlet boundary 
condition if)(0) = is relevant. For example, the half-line matrix Schrodinger equation 
(1.1) describes the behavior of n connected very thin quantum wires forming a one- vertex 
graph with open ends, and the boundary condition (1.6)-(1.8) imposes certain restrictions 
at the vertex. Such a problem is useful in designing elementary gates in quantum comput- 
ing and nanotubes for microscopic electronic devices, where, for example, strings of atoms 
may form a star-shaped graph. For details we refer the reader to [13,14] and the references 
therein. 

Our paper is organized as follows. In Section 2 we introduce various n x n matrix 
solutions to (1.1) that are needed later on. In Section 3 we introduce the Jost matrix J(k) 
and the scattering matrix S(k) and state some of their properties relevant to our study. 
In Section 4 we introduce two key transformations on the boundary condition (1.6)-(1.8) 
and analyze how those transformations affect the Jost matrix and the scattering matrix. 
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In Section 5 we provide the relevant properties of Jo(k) and So(k), the Jost and scattering 
matrices corresponding to V = 0; those properties are crucial in understanding similar 
properties when the potential in nonzero. In Section 6 we analyze the behavior of the 
Jost matrix at k = 0, which is needed in establishing Levinson's theorem. In Section 7 
we analyze the large-/c asymptotics of the Jost and scattering matrices. In Section 8 we 
analyze the bound states and the properties of the Jost matrix related to bound states. 
Finally, in Section 9 we establish Levinson's theorem, showing how the change in the phase 
of the determinant of the scattering matrix is related to the number of bound states. 

2. SCATTERING SOLUTIONS 

In this section we introduce certain n x n matrix solutions to (1.1) and recall some 
of their properties relevant to our study. We use C + to denote the upper-half complex 
plane and R for the real axis, and we let C+ := C + U R. Recall that we assume that the 
potential V appearing in (1.1) satisfies (1.3) and (1.5). 

The Jost solution to (1.1) is the n x n matrix solution satisfying, for k G C+ \ {0}, 
the asymptotics 



f(k, x) = e* kx [I n + o{l/x)}, f'{k, x) = ik e* kx [I n + o{l/x)}, x -> +oo, (2.1) 



and it is known [1,2] that f(k,x) and f'(k,x) are analytic in k G C + and continuous in 
k G C+ for each fixed x. We remark that f(k,x) corresponds to the quantity E(—k,x) 
described on p. 28 of [1]. From (2.1) it is seen that each of the n columns of f(k,x) 



The matrix Schrodinger equation (1.1) has the nxn matrix solution g(k, x) satisfying, 



where I n denotes the nxn identity matrix. It satisfies the integral equation 




(2.2) 



exponentially decays to zero as x — > +oo for each fixed k G C + . 
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for each k E C+ \ {0}, the asymptotics 

g(k,x) = e- tkx [I n + o(l/x)l g'(k,x) = -tke- tkx [I n + o(l/x)l x -> +00, (2.3) 

and g(k,x) corresponds to the quantity E^\— k, x) described on p. 28 of [1]. It is known 
[1] that g(k, x) and g'(k, x) are analytic in k E C + and continuous in k E C+ \ {0} for each 
fixed x. From (2.3) it is seen that each of the n columns of g(k 7 x) exponentially grows as 
x — > +00 for each fixed k E C + . 

As indicated on p. 28 of [1], for each k E C+ \ {0}, the combined 2n columns of 
f(k,x) and g(k,x) form a fundamental set of solutions to (1.1), and hence any column- 
vector solution u(k,x) to (1.1) can be written as 

u(k,x) = f(k,x)£ + g(k,x)ri, (2.4) 

for some constant column vectors £ and r\ in C n . 

The regular solution <p(k, x) is the nx n matrix solution to (1.1) satisfying the initial 
conditions 

(p(k,0) = A, <p'(k,0) = B, (2.5) 

where A and B are the matrices appearing in (1.6)-(1.8). For each fixed x E R + , it is 
known that <p(k, x) is entire in k in the complex plane C. Note that 

ip(-k,x) = (p(k,x), k E C, x E R + , (2.6) 

because k appears as k 2 in (1.1) and the initial values given in (2.5) are independent of k. 

We will use [F; G] := FG' — F'G for the Wronskian and use an asterisk to denote 
complex conjugation. It can be verified directly that for any two n x n solutions (f)(k, x) 
and ip(k,x) to (1.1), each of the Wronskians [<f>(k*,x)^;ip(k,x)] and [<f>(—k*,x)^;ip(k,x)] 
is independent of x. By evaluating the values of the Wronskians at two different x-values, 
say x = and x = +00, we can obtain various useful identities. For example, we have 

[f(±k,x)^j(±k,x)] = ±2ikl n7 kER, (2.7) 
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[f(-k*,xY;f(k,x)] = 0, keC+. (2.8) 

Let us add that, for each fixed x G R + , if a solution (f)(k, x) is analytic in k G C + then 
(f)(—k*, x)t becomes an analytic function of k in C + ; on the other hand, <f>(k*, x)^ becomes 
an analytic function of k in the lower- half complex plane C~. 

3. THE JOST MATRIX AND THE SCATTERING MATRIX 

In this section we introduce the Jost matrix J{k) and the scattering matrix S(k) for 
(1.1) with a selfadjoint matrix potential V in L\(R + ) and with the selfadjoint boundary 
condition (1.6)-(1.8). We also recall or establish some of their properties relevant to our 
study. 

The Jost matrix J(k) is defined in terms of a Wronskian as 

J(k) := [f(-k* , x)t; <p(k, x)], kE C+, (3.1) 

where f(k,x) is the Jost solution appearing in (2.1) and cp(k,x) is the regular solution 
appearing in (2.5). Since the Wronskian in (3.1) is independent of x, by evaluating its 
value at x = 0, with the help of (2.5) we get 

J(k) = f(-k\ 0) t S - f'(-k\ 0)U, (3.2) 

where A and B are the matrices appearing in (1.6)-(1.8). Note that the domain of J(k) is 
C+ because /(—&*, 0)"*" and f'(—k*,0y are analytic in k G C + and continuous in k G C+, 
as f{k, 0) and f'(k, 0) are analytic in k G C + and continuous in k G C+. 

We quote the following fundamental results from [2] regarding the Jost matrix and its 
inverse and their small- A; behavior. We refer the reader to Theorem 4.1 and Theorem 6.3 
of [2] for further details. These results are later needed in the derivation of Levinson's 
theorem. 

Theorem 3.1 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Then: 
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(a) J(k) is analytic in C + and continuous in C+. 



(b) As k — > in C+ we have 

J(k) = sp,- 1 



kAi + o(k) kB x + o(k) 
kd + o{k) V + o(l) 



(3.3) 



where P\ and Pi are some permutation matrices, and S is an invertible matrix, and 
Ai, Bi, C\ , and Vq are some constant matrices of sizes fix fx, fxx (n — /i), (n — fx) x fx, 
and (n — fx) x (n — fx), respectively, in such a way that Ai and Vq are both invertible. 
Here, fx is the geometric multiplicity of the zero eigenvalue of the zero-energy Jost 
matrix J(0). 

(c) J(k) is invertible for k e R \ {0}, and in fact J{k)~ x is continuous for k e R \ {0}. 



(d) As k — > in C+, we have 

r 1 

J(k)' 1 =SP 1 



k 



Ai\l» + o(l)] -Ai 1 B 1 Vq 1 + o(1) 



P2S- 1 . 



(3.4) 



l-V^dA^ + oil) V- x +o{l) 
Hence, J(k)~ x is either continuous at k = 0, or it has a simple pole at k = 0. In 
particular, kJ(k) -1 has a well-defined limit at k = 0. 

The following result is already known, but we provide a brief proof for the reader's 
benefit, as the information contained in the proof is relevant to our study. 

Proposition 3.2 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Then, 
the regular solution ip(k,x) can be expressed in terms of the Jost solution f{k,x) and the 
Jost matrix J(k) as 

V(K x) = ^ f(k, x) J(-k) - ^ f(-k, x) J(k), keK\ {0}. (3.5) 

PROOF: With the help of (2.1) we see that the combined 2n columns of f(k,x) and 
f(—k,x) form a fundamental set of column-vector solutions to (1.1) for any k E R \ {0}, 
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and hence we can write the regular solution <p(k, x) as 

<p(k, x) = f(k, x) dik) + f(-k, x) C 2 (k), keH\ {0}, (3.6) 

for some n x n matrices C\{k) and C2(k) depending only on k but not on x. Using (3.6) 
in the Wronskians [/(±fc, (p(k, x)], with the help of (2.7), (2.8), and (3.1) we obtain 

Ci(k) = ^J(~k), C 2 (k) = ~J(k) 1 

yielding (3.5) for real nonzero values of k. | 

Let us define the nx n physical solution ^(k,x) to (1.1) as 

V(k,x) := -2ikip(k,x) J(fc) -1 , (3.7) 

where <p(k, x) is the regular solution appearing in (2.5) and J(k) is the Jost matrix defined 
in (3.1). The scattering matrix S(k) is defined as [2,10-12] 

S(k):=-J(-k)J(k)-\ keK. (3.8) 

As we elaborate in Section 8, J(k) can uniquely be defined only up to a multiplication on 
the right by an invertible constant matrix. On the other hand, as seen from (3.8), such 
a postmultiplication of J(k) by an invertible matrix does not affect S(k). Hence, S(k) is 
uniquely determined by the potential V and the boundary condition (1.6), independently 
of the particular parametrization used in (1.6)-(1.8). In general, S(k) is defined only for 
real k because J(—k) in general cannot be extended from k G R to k G C + . The continuity 
of S(k) at k = has recently been established [2]. Even when J{k)~ x may not exist at 
k = 0, it has been shown [2] that the product on the right-hand side in (3.8) has a well- 
defined limit as k — )■ in R, and hence the domain of S(k) is k G R. The small-/c behavior 
of S(k) is quoted from [2] in the following and hence a proof is omitted. 

Proposition 3.3 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Then, 
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the scattering matrix S(k) defined in (3.8) is continuous for k G R including k = 0, and 
we have S(k) = S(0) + o(l) as k — > in R wif/i 







S(0) = SP 2 



-1 



P 2 S 



-1 



(3.9) 



-l 



where fi is the geometric multiplicity of the zero eigenvalue of the zero- energy Jost matrix 
J(0), P2 is an n x n permutation matrix, S is an n x n constant invertible matrix, Ai is 
a fx x fx constant invertible matrix, and C\ is an (n — fj) x 11 constant matrix. 

We note that the quantities fx, P2, S, Ai, and C\ appearing in (3.9) are the same as 
those appearing in (3.3) and (3.4). 

Proposition 3.4 Consider the self adjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Then 
for each x G R + the physical solution given in (3.7) is continuous for k G R and can be 
written as 



where S(k) is the scattering matrix defined in (3.8). 

PROOF: Using (3.5) and (3.8) in (3.7) we get (3.10) for k G R\ {0}. From the continuity 
of f(k,x) and S(k) for k G R including k = [2], it follows that (3.10) also holds at 
k = 0, and hence ty(k,x) is continuous in k G R for each fixed x G R + . We can verify 
the continuity of ^(k, x) in k G R in an alternate way. As stated in Theorem 3.1(d), even 
though J(fc) -1 may not exist at k = 0, the quantity kJ{k)~ l is continuous for k G R. We 
recall that (p(k, x) is entire in k for each x G R + . Thus, for each fixed x G R + , the physical 
solution ty(k,x) defined in (3.7) is continuous in k G R. | 

Some useful properties of the scattering matrix S(k) are provided in the following 
proposition. Thanks to the recent result [2] on the small-fc limit of S(k), the properties 
listed below hold for any real k, including k = 0. 



:=f(-k,x) + f(k,x)S(k), 



k G R, 



(3.10) 
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Proposition 3.5 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Then 
the scattering matrix S(k) defined in (3.8) is unitary for k G R and satisfies 

S(-k) = S(k)- 1 = S(fc)t, keK. (3.11) 

PROOF: Using (3.10) in the Wronskian x)t ; x)], with the help of (2.7) and (2.8) 
we get 

[tf^x^^x)] = -2ikl n + 2ikS(h)*S(k), keK. (3.12) 

On the other hand, using (3.7) in the same Wronskian, with the help of (1.7) and (2.5) we 
obtain 

[tf (fc, x) f ; (Jfe, x)] = -(2zA;) 2 [J(A;) t ]- 1 - fiU) J(fc)" 1 = 0, fc G R \ {0}. (3.13) 

In fact, (3.13) holds also at k = by letting k — > and noting that /cJ(/c) _1 has a well- 
defined limit [2] as fe — >■ in C+, as stated in Theorem 3.1(d). Comparing (3.12) and (3.13) 
we then get S(k)^S(k) = I n for k G R, which yields S(k)- 1 = S(fc) f for k G R. To establish 
S(—k) = S(ky for G R we proceed by evaluating the Wronskian A;, x)t; x)] in 
two different ways. First, using (3.10) in that Wronskian, with the help of (2.7) and (2.8) 
we get 

x)*; x)] = 2ik (S(k) - S(—k)^) , fc G R. (3.14) 

On the other hand, using (3.7) in the same Wronskian, with the help of (1.7), (2.5), and 
(2.6) we obtain 

[^(-k,x)^;^(k,x)} = {2ik) 2 [J{-k)^}- 1 (A^B-B^A) J(fc)" 1 =0, k G R\{0}. 

(3.15) 

For the same reason (3.13) holds at k = 0, we conclude that (3.15) also holds at k = 0. 
Thus, comparing (3.14) and (3.15) we conclude that S(-k) = S(fc)t for k G R. | 
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4. TRANSFORMATIONS 



In this section we remark how the Jost solution and the regular solution to (1.1), 
the Jost matrix, and the scattering matrix change if the matrices A and B used in the 
parametrization of the boundary condition (1.6)-(1.8) undergo a transformation without 
affecting (1.7) and (1.8). In particular, we consider a multiplication on the right by an 
invertible matrix, a unitary transformation by a unitary matrix, and a combination of 
those two transformations. The results will be useful in analyzing the large-/c asymptotics 
of various quantities and in the derivation of Levinson's theorem. 

Proposition 4.1 Consider the self adjoint matrix Schrddinger operator with the self adjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Let 
A and B be the matrices appearing in (1.6)-(1.8), f(k,x) be the Jost solution to (1.1) 
satisfying (2.1), <p(k,x) be the regular solution to (1.1) satisfying (2.5), J{k) be the Jost 
matrix defined in (3.1), and S(k) be the scattering matrix defined in (3.8). Then: 

(a) Under the transformation V (->■ V and (A,B) >->■ (AT, ST), where T is an invertible 
n x n matrix, we have 

(f,<p,J,S)^(f,<pT, JT,S). 

(b) Under the unitary transformation V H> MWM and (A,B) H> {M^AM.M^BM), 
where M is a unitary n x n matrix, we have 

(f, ip, J, S) h-> (M^fM, MVM, M f JM, M^SM). 

(c) Under the unitary transformation V >->■ M^VM with a unitary matrix M and the com- 
bination of three consecutive transformation (A, B) h-> (M^ ATiMT 2 , BTiMT 2 ) 1 
first by a right multiplication by an invertible matrix Ti, then by the unitary trans- 
formation with M, followed by a right multiplication by an invertible matrix T 2 , we 
have 

(f, ip, J, S) ^ (Mt/M, M 1r (pT 1 MT 2 , M^JT 1 MT 2 , M ] SM). (4.1) 
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PROOF: The proof is obtained by direct verification and by checking that the boundary 
condition (1.6)-(1.8), the Schrodinger equation (1.1), and the relevant conditions and def- 
initions in (2.1), (2.5), (3.1), and (3.8) all remain satisfied. Finally, the transformation in 
(c) is obtained from the results in (a) and (b). | 

We note that the transformation V i— > V and (A, B) h-> {AT, BT) with an invertible 
matrix T is just a change of parametrization in the boundary condition (1.6)-(1.8). On 
the other hand, the unitary transformation V H> M^VM and (A, B) ^ (M^AM, M^BM) 
with a unitary matrix M is a change of representation in the sense of quantum mechanics. 

Motivated by the general selfadjoint boundary condition [3,19,20] in the scalar case, 
i.e. the case with n = 1, we are interested in going from the pair A and B appearing in 
the selfadjoint boundary condition (1.6)-(1.8) to the special pair A and B, where we have 
defined 

A := — diag{sin^!, . . . ,sin# n }, B := diag{cos#i, . . . ,cos# n }, (4.2) 

with the real parameters 9j taking values in the interval (0, 7r]. The special case 9j = tv 
corresponds to the Dirichlet boundary condition and the case 9j = n/2 corresponds to the 
Neumann boundary condition. We assume that there are tin values with 9j = n/2 and no 
values with 9j = n, and hence there are hm remaining values, with um '■= n — — rip, 
such that those ^-values lie in the interval (0, ir/2) or (n/2,n). Our analysis takes into 
consideration the special cases where any of um, tid, and Nm may be zero or n. In our 
notation the subscripts M, D, and refer to "mixed," "Dirichlet," and "Neumann," 
respectively. We assume that the 6>j-values in (4.2) are ordered in such a way that the first 
um values of 9j correspond to the mixed conditions, the next n d values correspond to the 
Dirichlet conditions, and the remaining tin values correspond to the Neumann conditions. 

We will provide the explicit steps to go from any pair of matrices A and B satisfying 
(1.6)-(1.8) to the pair A and B given in (4.2) and yet still satisfying (1.6)-(1.8) with A and 
B replacing A and £?, respectively, there. For this, we need some auxiliary results. 
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Starting with A and B satisfying (1.6)-(1.8), let us define 

E := (A^A + B^B) 1 / 2 , (4.3) 
so that E is positive, and hence E is uniquely defined. 

Proposition 4.2 Let A and B be a pair of matrices satisfying (1.6) and (1.7), and let E 
be the matrix defined in (4-3). Then: 

(a) The matrix E is invertible and satisfies 

E = E\ E- 1 (A t A + B j B)E- 1 =I n . (4.4) 

(b) We have 

(B ± iA)E~ 2 (B Jf =f iA^) = I n , (4.5) 
and hence the matrices (B ± iA) and (B^ ± iA^) are all invertible and in fact 

(B ± iA) -1 = E~ 2 (B Jf =f iA^). (4.6) 

(c) The matrix U defined as 

U := (B — iA)E~ 2 (B^ — iA^), (4.7) 

is unitary, and hence it satisfies UU> = UW = I n . 

(d) The matrix U defined in (4-7) can also be written as 

U = (B -iA)(B + iA)~\ (4.8) 

and hence from = U~ x it follows that 

= (B + iA)(B -iA)' 1 . 
15 



PROOF: The proof of (a) readily follows from (4.3). To prove (b) we let 

'BE' 1 AE' 1 

C:= 

AE' 1 -BE' 1 

and, by using (1.7) and (4.4), we directly verify that C^C = I 2n , proving the unitarity of 
C. We must then also have CC^ = ^n, implying 

AE~ 2 A ] + BE~ 2 B^ = 7 n , BE~ 2 A ] - AE~ 2 B^ = 0. (4.9) 

With the help of (4.9), from the identity 

(B ± iA)E- 2 (B* =f zA f ) = (AE~ 2 A^ + BE~ 2 B^) t i(BE~ 2 A^ - AE~ 2 B^), 

we obtain (4.5) and hence (b) is proved. Let us now turn to the proof of (c). With the help 
of (4.3) and (4.5) we directly verify that the matrix U defined in (4.7) satisfies UU^ = I n , 
and hence (c) is proved. From (4.6) and (4.7) we get (4.8) and hence (d) is also proved. | 

Proposition 4.3 Let A and B be a pair of matrices satisfying (1.7) and (1.8), and let A 
and B be the matrix pair appearing in (4-2). We then have 

A = M^AT 1 MT 2l B = M ] BT 1 MT 2 , (4.10) 

for some unitary matrix M and for some invertible matrices T\ and T 2 . 

PROOF: We can diagonalize the unitary matrix U appearing in (4.7) and (4.8) by using 
a unitary matrix M so that 

M^UM = diag{e 2lCl , . . . , e 2lC ™ }, (4.11) 

where the constant parameters satisfy Q G (0, n]. Let us define 

Y :=diag{e lCl ,...,e lC "}. (4.12) 

With the help of a permutation matrix P, we can reorder Q as 9j in the manner described 
below (4.2), namely, the first um values of 9j lie in (0, ir/2) U (n/2, 7r), the next no values 
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of 9j are all equal to tt, and the remaining values of 9j are all equal to tv/2. Thus, from 
(4.11) and (4.12) we obtain 

M^UMP = Y 2 P = diag{e 2 ^ , . . . , e 2i9 ~ }, 

yP = diag{e l V..,e i9 ™}, (4.13) 
Y~ X P = diagje"^ 1 , . . . , e' 10 - }. 
On the other hand, from (4.2) we see that 

S-zI = diag{e^ 1 ,...,e^™}, (4.14) 

and hence 

(B-iA)- 1 = S + zI = diag{e- ie \...,e-^"}. (4.15) 

We are now ready to prove (4.10). The invertibility of (B + iA) is assured by (4.6) and 
hence we have 

(B + iA)(B + iA)' 1 =I n . (4.16) 

From (4.8) we have 

(B - iA)(B + iA)' 1 = U. (4.17) 

Let us premultiply (4.16) and (4.17) by and postmultiply them by MY~ l P in order 
to obtain 

M\B + iA)(B + iA)- x MY- x P = Y~ X P, (4.18) 

M\B - iA)(B + iA)- 1 MY~ 1 P = M^UMY~ X P, (4.19) 

where we have used the unitarity property M^M = I n . From (4.11) and (4.12) we know 
that MWM = Y 2 and hence M^UMY~ l P = YP. Thus, we can rewrite (4.19) as 

M j (B-iA)(B + iA)- 1 MY- 1 P = YP. (4.20) 

On the other hand, from (4. 13)- (4. 15) we see that 

YP = B — iA, Y~ X P = B + iA. (4.21) 
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Using (4.21) we can rewrite (4.18) and (4.20) as 

M\B + iA)(B + iA)- l M{B + iA) = B + iA, (4.22) 

M\B - iA)(B + iA) _1 M(S + iA) = B- iA. (4.23) 

Letting 

Ti := (B + iA)-\ T 2 :=B + iA, (4.24) 
because of (4.6) and (4.15) we observe that Ti and T 2 are invertible and in fact 

T± 1 = B + iA, T 2 ~ l =B-iA. (4.25) 

Using (4.24) we can rewrite (4.22) and (4.23), respectively, as 

M f (B + iA)T 1 MT 2 = B + iA, (4.26) 

M\B — iA)T 1 MT 2 = B — iA. (4.27) 
By subtracting and adding, respectively, from (4.26) and (4.27) we obtain (4.10). | 

Because of Propositions 4.1 and 4.3, there is not much loss of generality in using the 
special boundary parametrization with A and B given in (4.2). The relevant results can 
then be transformed to obtain the corresponding results in the parametrization with A 
and B appearing in (1.6)-(1.8). Let us use a tilde to denote the quantities obtained under 
the transformation given in Proposition 4.1(c). With the help of (4.1) and (4.10), we 
can obtain the corresponding quantities in any boundary parametrization with A and B 
satisfying (1.7) and (1.8). In other words, when we have 

V(x) =MV(x) M\ A = MAT2 1 M t Tj~ 1 , B = MBT 2 ~ l M ] T^ 1 , (4.28) 



we then get 

f(k,x) = Mf{k,x)M\ J(k) = MJ(k)T 2 ~ 1 M Ji T^ 1 , S(k) = MS{k)M\ (4.29) 
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where M is the unitary matrix appearing in (4.11) and T x 1 and T 2 1 are the matrices 
specified in (4.25). 

5. THE JOST AND SCATTERING MATRICES WITH ZERO POTENTIAL 

In order to understand the large-/c behavior of the Jost matrix J(k), its inverse J(/c) _1 , 
and the scattering matrix S(k), we need to understand those behaviors when the potential 
V appearing in (1.1) is identically zero. In that case, let us use the subscript to denote the 
corresponding quantities and write Jo(k) and So(k) for the Jost and scattering matrices, 
respectively, corresponding to V = 0. 

When V = 0, from (2.2) we get f(k, x) = e lkx I n , and hence (3.2) and (3.8) yield 

J (k)=B-ikA, Mk)- 1 = (B-ikA)-\ S (k) = -(B + ikA)(B -ikA)' 1 . (5.1) 

In the boundary parametrization with A and B in (4.2), the corresponding quantities are 
given by diagonal matrices, where we have 

Jo(k) = B - ikA = diag{cos#i + ik sin#i, . . . , cos6 nM + z/csin# nM , —I nD ,ikI nN }, (5.2) 

Jo(k) diag | cos( ^ _|_ ^ g ^ n ^ 1 ' • •• ' cos Q nM _|_ ik sin 9 nM ' ^ nrn ik ^ nN } ' 

~ f -cos^i + z/csin^i - cos 9 nM +iksm6 nM \ 

S {k) = diag <^ • • • , ~ a ,.,•/) » -In D Jn N \ ■ (5.4) 

L cosyi + z/csmPi cosB nM + iksmv nM J 

Then, from (5.3) and (5.4) we see that, as — > oo in C, 

(k) _1 = diag{0 nM , -I nD , niV } + — diag{csc X , . . . , csc 6> nM , nD , I njv } + 0(l//c 2 ) , (5.5) 

So(k) = Z + jZ 1 +0(l/k 2 ), (5.6) 

where we have defined 

Z := diag{/ nM , -I nD ,I nN }, Z x := diag{cot6»i, . . . ,cot0 nM ,O n£) ,O njv }, (5.7) 
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with Oj denoting the j x j zero matrix. 

The results in the following proposition are needed in Section 7. 

Proposition 5.1 The Jost matrix Jo(k) appearing in (5.1) is invertible when k — > oo in 
C. The matrix Jo(/c) _1 and the scattering matrix So(k) given in (5.1) satisfy, as k — >■ oo 
in C, 

Jo(fc)- 1 =T 1 MT 2 diag{0 nM ,-I nD ,0 nN }M^ + 0(l/k), (5.8) 

AJo(k)- 1 = ^Mdiag{-I nM ,O nD ,-I nN }M^ + 0(l/k 2 ), (5.9) 

So(k) = S (oo) + MZiM^ + 0(l//c 2 ), S (oo) := MZ M^, (5.10) 

where T\ and T 2 are the matrices in (4-24), M is the unitary matrix in (4-11), A and B 
are the matrices appearing in (1.6)-(1.8), A and B are the matrices defined in (4-2), and 
Zq and Z\ are the matrices defined in (5.7). Thus, as k — > oo in C we have 

Jo(k)- 1 = 0(1), AJo(k)- 1 = 0(l/k), S (k) = So(oo) + 0(l/k). (5.11) 

PROOF: Exploiting the unitarity properties M _1 = M f and M f M = I n , from (4.28) and 
(4.29) when V = we get 

J^k)- 1 =T 1 MT 2 J Q {k)- 1 M\ AJoik)- 1 = MAJ Q {k)- 1 M\ S (k) = MS (k)M*. 

(5.12) 

From (5.3) we conclude that Jo(k) is invertible when k — > oo in C, and hence the first 
equality in (5.12) implies that Jo(/c) _1 exists when k — > oo in C. Using (4.2), (5.5), and 
(5.6) in (5.12), we get the expansions (5.8)-(5.10) as k — > oo in C. | 

6. SMALL-A; BEHAVIOR 

The analysis of (1.1) near and at k = deserves a separate attention and we refer the 
reader to [2] for such an analysis. Under the assumption that the potential V satisfies (1.3) 
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and (1.5), we recall certain useful properties that will be needed in establishing Levinson's 
theorem. 

When k = from (1.1) we obtain the zero-energy matrix Schrodinger equation 

if;" = V(x) iGR + . (6.1) 

It is known [1,2] that the Jost solution f(k,x) to (1.1) appearing in (2.1) satisfies (6.1) if 
we replace k by in f(k,x), and that f(0,x) satisfies [1,2] 

/(0,x) = / n + o(l), / / (0,x) = o(l/x), x^+oo. (6.2) 

It is also known [1,2] that g(0,x), obtained by replacing k with in the matrix- valued 
function g(k,x) appearing in (2.3), is a solution to (6.1) and it satisfies 

g(0,x) = x[I n + o(l)}, g'(0,x) = I n + o(l), x -> +oo. (6.3) 

Thus, from (6.2) and (6.3) we see that the combined 2n columns of /(0, x) and g(0, x) form 
a fundamental set of solutions to (6.1), and hence (2.4) is valid even when k = 0. From 
(6.2) we see that the n columns of f(0,x) form n linearly independent solutions to (6.1) 
that remain bounded as x — > +oo. Similarly, (6.3) indicates that the n columns of g(0, x) 
form n linearly independent solutions to (6.1) that become unbounded as x — > +oo. 

On the left-hand side of (2.4) with k = 0, let us use a linear combination of n columns 
of the zero-energy regular solution cp(0, x), which itself is an n x n matrix solution to (6.1), 
and let us express such a column- vector solution as linear combinations of the 2n combined 
columns of /(0, x) and g(0,x), i.e. we let 

<p(0,x)u = f(0,x)Z + g(0,x)ri, (6.4) 

for some nonzero constant column vector u G C n . We know that the left-hand side in (6.4) 
satisfies (1.6) for any u G C n because <p(0,x) itself satisfies (1.6). We are interested in 
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knowing how many linearly independent bounded column-vector solutions to (6.1) we can 
form by using linear combinations of n columns of <p(0, x). Equivalently, we are interested 
knowing how many of the n linearly independent bounded column-vector solutions to (6.1) 
also satisfy (1.6). 

There are two possibilities for (6.4) with a nonzero column vector u G C n ; either 
ip(0, x) u is bounded as x — > +oo, in which case we must have £ 7^ and rj = 0, or cp(0, x) u 
is unbounded as x — > +00, in which case we must have rj ^ 0. In fact, the related results 
are known and quoted from [2] in the following proposition. 

Proposition 6.1 Consider the self adjoint matrix Schrddinger operator with the self adjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Then, 
we have 

(a) The nonzero column vector u G C n is an eigenvector of the zero-energy Jost matrix 
J(0) with the zero eigenvalue, i.e. u G Ker[J(0)] if and only if <p(0, x) u is bounded 
for x G R + . 

(b) For any column vector u in Ker[J(0)] there exists a unique column vector I; in Ker[J(0)^] 
such that 

<p(0,x)u = f(0,x)£. 
The map from Ker[J(0)] to Ker[J(0y] is a bijection. 

As in Theorem 3.1 let us use \i to denote the geometric multiplicity of the zero eigen- 
value of J(0). From Proposition 6.1 it follows that we can form exactly ji columns by 
using linear combinations of n columns of the zero-energy regular solution <p(0, x) in such 
a way that those fx columns form linearly independent solutions to (6.1) and they remain 
bounded as x — > +00. Furthermore, each of such \i column-vector solutions to (6.1) can 
also be expressed as a linear combinations of columns of /(0, x). In that sense, the integer \i 
indicates the maximal number of linearly independent bounded solutions to (6.1) that also 
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satisfy (1.6), and hence \x acts as a "degree" of the exceptional case [2] for the Schrodinger 
equation (1.1) with the boundary condition (1.6)-(1.8). In the purely generic case, i.e. 
when /j, = 0, the 2n combined columns of <p(0,x) and of f(0,x) are all linearly indepen- 
dent. In that case, each column of cp(0, x) can be expressed as a linear combination of n 
linearly independent columns of g(0,x). In the purely exceptional case, i.e. when \i = n, 
each column of <p(0,x) can be expressed as a linear combination of n linearly independent 
columns of /(0, x). 

From (3.3), we have the following conclusion. 

Corollary 6.2 Consider the selfadjoint matrix Schrodinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Then, 
the determinant of the Jost matrix J(k) defined in (3.1) has the small-k behavior 

det J(k) = ci/c M [l + o(l)], fc^0inC+, (6.5) 

where fx is the geometric multiplicity of the zero eigenvalue of the zero- energy Jost matrix 
J(0) and c\ is a nonzero constant. In fact, the value of c\ is given by 

ci := (det Pi) (det P 2 ) (detail) (det I?o), 

where P\ and are the n x n permutation matrices appearing in (3.3) and hence their 
determinants are either 1 or —1, and A\ and Vq are the invertible matrices appearing in 
(3.3) and hence their determinants are nonzero. 

From (3.9) we see that the zero-energy scattering matrix S(0) has only two eigenvalues, 
namely +1 and —1. In the next proposition we prove that the eigenvalue +1 of S(0) has 
multiplicity (both geometric and algebraic) equal to fi and that the eigenvalue —1 has 
multiplicity (both geometric and algebraic) equal to n — fx. Thus, the value of fi is uniquely 
determined from S(0) alone. 

Proposition 6.3 Consider the selfadjoint matrix Schrodinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Let J(k) 
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and S(k) be the corresponding Jost matrix and the scattering matrix defined in (3.1) and 
(3.8), respectively. Then: 

(a) S(0) has two eigenvalues, which are +1 and — 1. 

(b) The geometric and algebraic multiplicities of the eigenvalue +1 of S(0) are both equal 
to f/,, which is the geometric multiplicity of the zero eigenvalue of J(0). 

(c) The geometric and algebraic multiplicities of the eigenvalue — 1 of S(0) are both equal 



PROOF: By (3.9), we see that S(0) is similar to a lower-triangular matrix whose diagonal 
entries coincide with the diagonal entries of the diagonal matrix diag{/ M , — J n _ M }. Since a 
similarity transformation does not change the eigenvalues, we conclude that fS'(O) has the 
eigenvalue +1 with the algebraic multiplicity fx and the eigenvalue —1 with the algebraic 
multiplicity n — fx. We thus only need to determine the geometric multiplicity for each 
eigenvalue. By (3.11), we know that S(0) is unitary and hence it can be diagonalized 
with a unitary matrix M\. Thus, S(0) = MiDMl for some diagonal matrix D, with \x of 
the diagonal entries being +1 and n — fx of them being —1. The algebraic and geometric 
multiplicities of eigenvalues remain invariant under a similarity transformation with a uni- 
tary matrix, and furthermore the algebraic and geometric multiplicities of each eigenvalue 
of a diagonal matrix are equal to each other. Thus, the geometric multiplicities of the 
eigenvalues +1 and —1 are given by fx and n — fx, respectively. | 

7. LARGE-*; BEHAVIOR 

In establishing the large-/c behavior of various quantities related to (1.1), it is sufficient 
for the potential to satisfy (1.2) rather than the stronger condition (1.5). The following 
matrices will be useful in our large- k analysis: 



to n — fx. 




(7.1) 
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-I poo l>Z -I i* OO l>Z 

Qs:= 4j dz J q d V V ( z ) V (y), Q4(k):=-j^ dz dye 2lhz V(z)V(y), (7.2) 

-| POO pz -i poo pz 

Q 5 (k):=-J dzj dye 2 ^V(z)V(y), Q 6 (k) := - J dz J dy e 2 ^V(z)V(y). 

(7.3) 

We emphasize that Q± and Q3 are independent of k while the remaining four matrices are 
functions of k. 

Proposition 7.1 Assume that V in (1.1) belongs to L 1 (R + ). Then, Q\ and Q 3 are well 
defined constant matrices. Furthermore, each of the four matrix quantities Q2(k), Q 4 (k), 



Qsik), and Qe(k) is well defined for k G C+, and each one of them has the behavior of 



o(l) as k — )> 00 in C+. 



PROOF: The integrals defining Q\ and Q3 exist because V G L 1 (R + ). For k G C+, the 
coefficients of 2ik in the exponents appearing in Q 2 (^) 5 Q^ik), Qs(k), and Qe(k) are all 
nonnegative, and hence each of those exponential terms are bounded by one in absolute 
value. Furthermore, V G L 1 (R + ), and hence we have the estimate 

PX px P 00 

/ dzV(z) < / ds|W*)||< / dz\\V(z)\\, 
Jo Jo Jo 

where the integral on the right-hand side converges. Thus, we can establish that the 



integrals used in defining Q2(k), Q 4 (k), Q*,(k), Qe(k) all exist for k G C+. Finally, the 



behavior of o(l) as k — > 00 in C+ for Q 2 (^) 5 Q±{k), Qs(k), Qe(k) is established with the 
help of the Riemann-Lebesgue lemma on the appropriate integrals. | 

Proposition 7.2 Consider the matrix Schrddinger equation (1.1) with the self adjoint po- 
tential V satisfying (1.2) and (1.3). Let f(k, x) be the corresponding Jost solution satisfying 



(2.1). Then, we have as k — >■ 00 in C+ 

f(-k*, 0) t = I n +^ [-Qi + Q 2 (k)]+^ [-Q 3 - Q 4 (k) + Q 5 (k) + Q 6 (k)]+0(l/k 3 ), (7.4) 

f'{-k*, 0)t = ikl n - Qi- Q 2 (k) + ^[Qs- Q 4 (k) + Q 5 (k) - Q 6 (k)] + 0(l/k 2 ), (7.5) 
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where Qi, Q2(k), Q3, Q^k), Q$(k), Q&(k) are the matrices defined in (7.1)-(7.3) with the 
properties outlined in Proposition 7.1. 

PROOF: Writing the Jost solution f(k,x) in terms of m(k,x) := e~ lkx f(k,x), from (2.2) 
we obtain the integral relations 

™ 00 



m(k,x) = I n + ^ / dy 



e 2ik{y-x) _ y 



V(y)m(k,y), 



m'(k,x) = - / dye 2ik ^- x) V{y)m(k,y). 

J X 

Iterating (7.6) and (7.7), we obtain 



(7.6) 
(7.7) 



1 [°° 
m(k, x) =I n + — / dy 



+ 



1 



(2iky 



dy 



e 2ik(y-x) _ ^ 
e 2ik(y-x) _ ^ 



V{y) 
V(y) f 

J V 



dz 



e 2ik{y-z) _ y 



V(z)+0(l/k 3 ), 
(7.8) 



m'(/c, x) = 



dy e 



2ik(y—x) 



V(y) 



1 



- ttt / d y e 



D 2ik(y-x) 



r°° r 1 
V{y) dz e 2lk{y - z) -1 \V(z) + 0(l/k 2 ). 



(7.9) 



We evaluate (7.8) and (7.9) at x = 0, and we rewrite the double integrals in them by 
changing the order of integration. Then, with the help of 



f(k, 0) = m(k, 0), f'(k, 0) = ik m(k, 0) + m'(k, 0), 



(7.10) 



we obtain the expansions for f(k, 0) and f'(k, 0) as k — > 00 in C+. Finally, by replacing k by 
—k*, taking the adjoints, and using [V(y) V(z)} ] = V(z) V(y), with the help of (7.8)-(7.10) 
we obtain (7.4) and (7.5). | 

Proposition 7.3 Consider the Schrddinger operator corresponding to (1.1) with the self- 
adjoint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.2) and (1.3). 
Let J(k) be the corresponding Jost matrix defined in (3.1). Then, we have as k — >■ 00 in 



C+ 



J(k) = -ikA + B + [Qx + Q 2 (k)]A + \p{k) + 0(l/k 2 ), 



(7.11) 
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where we have defined 



P(k) := [-Qx + Q 2 (k)]B + [-Q 3 + Q A {k) - Q 5 {k) + Q 6 (k)]A, 

with A and B being the matrices appearing in (1.6)-(1.8), and Qi, Q2(k), Q3, Q±{k), 
Qb{k)i Qe(k) being the matrices in (7.1)-(7.3). 

PROOF: Using (7.4) and (7.5) in (3.2) we obtain (7.11). | 

Proposition 7.4 Consider the Schrddinger operator corresponding to (1.1) with the self- 
adjoint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.2) and (1.3). 
Let J(k) and Jo(k) be the corresponding Jost matrices appearing in (3.1) and (5.1), respec- 
tively, Qi and Q 2 (k) be the quantities in (7.1), and So (00) be the constant matrix defined 
in (5.10). Then, as k — > 00 in C+ we have 

J(k) Jo(k)- 1 = I n - 1 [Q x + Q 2 (k) S (oo)} + 0(l/k 2 ), (7.12) 
Jo(fc) J(k)- 1 =l n + ±-[Q 1 + Q 2 (k) So(oo)] + 0(l/k 2 ). (7.13) 

PROOF: By replacing B — ikA by Jo(k), as given in (5.1), and replacing B + ikA by 
Jo(— k), from (7.11) we get, as k — > 00 in C+, 

J(k) = Mk)~ QiJo(k)+^ Q 2 (k) J (-k)+^[-Q 3 +Q 4 (k)-Q 5 (k)+Q 6 (k)}A+O(l/k 2 ). 

(7.14) 

The invertibility of Jo(k) as k — > 00 in C+ is assured by Proposition 5.1. Let us multiply 
(7.14) on the right by J (/c) _1 and use J (—k) J (/c) _1 = —S (k), which follows from (3.8). 
We then obtain, as k — > 00 in C+, 

J(k) Jo(k)- 1 =/ n - 1 Q x - 1 Q 2 (k) S (k) 

+ \[-Qs + Q4k) - Qa(k) + Q 6 (k)]A J (/c) _1 + 0(l/k 2 ) J (/c) _1 . 

(7.15) 
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Using (5.11) in (7.15), we get 

J(k) Jo(k)- 1 = / n - 1 - 1 Q 2 (k) S (k) + 0(l/k 2 ), k -> oo in C+, 

which is also equivalent to (7.12) because of the third estimate in (5.11). The expansion 
in (7.13) is obtained by changing the signs of the 0(l//c)-terms in (7.12). | 

Proposition 7.5 Consider the Schrddinger operator corresponding to (1.1) with the self- 
adjoint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.2) and (1.3). 
Then, the corresponding Jost matrix J(k) defined in (3.1) satisfies 



J(k) = J (k)[I n + O(l/k)], £;^ooinC+, (7.16) 



det J(k) =c 2 k nM+nN [l + 0(l/k)], £;^ooinC+, (7.17) 

where Jo(k) is the matrix in (5.1), c 2 is a nonzero constant, and hm and un are the 
nonnegative integers defined after (4-2) and appearing in (5.2)-(5.4). 

PROOF: Note that (7.16) is apparent from (7.15). With the help of the first equality in 
(5.12) we get 

Jo(fc) = MJ (k) T^A^Tf 1 , (7.18) 

where T\ and T 2 are the invertible matrices appearing in (4.24) and (4.25), M is the unitary 
matrix appearing in (4.11), and Jo(k) is the matrix in (5.2). Using (5.2) and (7.18) we 
obtain 



det Jo(fc) = c 2 k nM+nN [l + 0(l/k)\, k ->■ oo in C+, (7.19) 
where we have defined 

C2:= Lil_4!L TTsin^, 

det[TiT 2 ] f\ 3 

Note that c 2 is well defined and nonzero because Ti and T 2 are invertible matrices and 
sin 9j 7^ for j = 1,...,um- The latter follows from the fact that those 9j all lie in 
(0,7r/2) U (tt/2,7t), as stated above (4.13). From (7.16) and (7.19) we then get (7.17). | 
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Next we present the large-/c asymptotics of the scattering matrix. 

Theorem 7.6 Consider the Schrddinger operator corresponding to (1.1) with the self ad- 
joint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.2) and (1.3). 
Then the corresponding scattering matrix S(k) defined in (3.8) satisfies: 

S(k) = S (oo) + ^ffl +0(l/k 2 ), k^±oo, (7.20) 

IK 

where G(k) is the matrix defined as 

G{k) := -2MZ 1 M^ + Q 1 S (oo) + S (oo) Q 1 + S (oo)Q 2 (k)S (oo) + Q 2 {-k), 

with M being the unitary matrix in (4- 11), Z\ the matrix in (5.7), So(oo) the matrix in 
(5.10), and Qi and Q 2 {k) the matrices in (7.1). 

PROOF: With the help of (3.8) we see that 

S(k) = J(-k) Joi-k^Soi^Joik) J(k)-\ (7.21) 

where Jo(k) and S (k) are the matrices defined in (5.1) and J(k) is the Jost matrix 
appearing in (3.1). Using (7.12) and (7.13) in (7.21), with the help of the identity 
So(—k)So(k) = I n , which follows from (3.11), we obtain 

S(k) = S (k) + —-^ + O(l/k 2 ), k^±oo, (7.22) 
%k 

where we have defined 

H(k) := QiSo(fc) + So(Wi + S (k)Q 2 (k)S (k) + Q 2 (-k). (7.23) 
Finally, using (5.10) in (7.22) and (7.23), we obtain (7.20). | 

8. BOUND STATES 

A bound state for the Schrddinger equation (1.1) with the boundary condition (1.6)- 
(1.8) at x = corresponds to a square-integrable column- vector solution satisfying (1.6). 
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Because of (1.3) and (1.6)-(1.8), the corresponding Schrodinger operator is selfadjoint, and 
hence a bound state must occur at a real value of k 2 . From (2.1), (2.3), and (2.4), we see 
that we cannot have any bound states at any positive values of k 2 because none of the 
combined 2n linearly independent columns of f(k,x) and g(k,x) can be square integrable 
on x E R + when k is real. Similarly, from (6.2) and (6.3) we see that none of the 
combined 2n linearly independent columns of f(0,x) and g(0,x) can be square integrable 
on x G R + , and since (2.4) also holds at k = 0, we can conclude that there cannot be a 
bound state when k = 0. Thus, a bound state, if it exists, can only occur when k 2 < 0, 
which corresponds to a value of k on the positive imaginary axis in C. 

Let us assume that k = in for some positive k, corresponds to a bound state, and let 
oj(iK,x) be a square-integrable column-vector solution satisfying (1.6). Then, (2.4) must 
hold at k = ik with ij = and for some nonzero column vector £ e C n , yielding 



u(ik, x) = f(in,x)£. 



(8.1) 



Let us show that £ must belong to the kernel of J(ik)^ because u(ik,x) must satisfy (1.6). 
Note that, from (1.6) and (8.1) we get 



Bt/(i«,0)£ + At/'(iK,0)£ = 0, 



(8.2) 



which is equivalent to 



[-<p'(iK, 0)t/(m, 0) + <p(i K , 0)t/'(i«, 0)] £ = 0, 



(8.3) 



or equivalently, in terms of the Wronskian, (8.3) can be written as 



[f(iK, xy,(p(iK,x)y£ = 0. 



(8.4) 



Comparing (3.1) and (8.4), we see that (8.4) is equivalent to 



J(ik)^ = 0, 



(8.5) 
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and hence £ belongs to Ker[J(m)t]. Thus, the determinant of J (in) must be zero. 

Conversely, consider any column vector of the form f(iK,x)£, where k = in corre- 
sponds to a zero of det J(k) on the positive imaginary axis and £ G C n is a nonzero 
column vector belonging to Ker[J(m)t]. Then, /(m,x)£ must be a bound-state column- 
vector solution to the corresponding Schrodinger operator. To verify this, we must prove 
that f(iK, x) £ is a solution to (1.1), it is square integrable on x G R + , and it satisfies the 
boundary condition (1.6)-(1.8). It is clearly a solution to (1.1) because f(k,x) is a n x n 
matrix solution to (1.1). It is square integrable because f(k,x) exponentially decays to 
zero for each k G C + as x — > +oo, as apparent from (2.1). Finally, it satisfies the boundary 
condition (1.6) because 

-Bif(i K , 0)£ + A^f'(iK, 0)£ = J(ik)^ = 0, 

as seen from (8.2)-(8.5). The multiplicity of the bound state at k = ik is equal to the 
dimension of the kernel of J{in)\ which is also equal to the dimension of the kernel of 
J{ik). 

Let us now show that a bound-state column-vector solution at k = in must have the 
form cp(iK, x)a for some constant nonzero column vector a G C n that belongs to the kernel 
of J {in) in such a way that 

<f(iK,x)a = f(iK,x)P, (8.6) 

where (3 G Ker[J(m)']. In other words, we must show that (p(in,x)a with a G Ker[J(m)] 
satisfies (1.6), is square integrable, and satisfies (8.6) for some column vector /3 G Ker[J(iK)" t "] 
Note that (1.6) is satisfied because with the help of (1.7) and (2.5) we get 

°)" + A*<p'(iK, 0)a = (-B ] A + A ] B)a = 0. 

Next, let us show that <p{iK,x)a exponentially decays to zero as x — > +oo, and hence it is 
square integrable on x G R + . From (2.4) we see that Lp(iK,x)a can be written as a linear 
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combination of the 2n linearly independent columns of f(m, x) and g(in, x), i.e. there exist 
some constant column vectors (3 and 7 in C n such that 

<p(iK, x) a = f(iK, x) (3 + g(iK,x) 7, x G R + . (8-7) 

Let us now evaluate the Wronskian-related quantity [f(in, x) 1 "; (p(in, x)]a using (8.7). With 
the help of (3.1) we get 

[f(in, x)^; ip(iK, x)]a = J(in)a = 0, (8-8) 

because a e Ker[J(w)]. On the other hand, using (8.7) we get 

[/(*«, a:) 1 ; a)]a = [/(*«, ^) f ; /(*«, x)/3 + g(iK, x) 7] 

= [/(*«, x) f ; x)] /? + x) f ; x)] 7 

(8.9) 

= [/(««, ^) f ;^(^,^)] 7 

= 2^7, 

where we have used (2.1) and (2.3) to evaluate the relevant Wronskians. Comparing (8.8) 
and (8.9) we see that 7 = and hence (8.6) is satisfied for some nonzero column vector /3 
in C n . Because of (2.1), from (8.6) we conclude that (p(in, x)a decays exponentially to zero 
as x — > +00 and hence it is square integrable. Note that, /3 must belong to Ker[J(m)t] as 
a result of our earlier argument that if f(in, x) f3 is a bound state then f3 must belong to 
the kernel of J(m)L 

Let us emphasize that (8.6) establishes a bijection a (->■ ft between Ker[J(m)] and 
Ker[J(m)T] for any k = in that is a zero of det J(k) on the positive imaginary axis. 
Since f(iK, x)£ corresponds to a bound state with £ e Ker[J(m)" 1 "], from (2.1) we conclude 
that there are as many linearly independent bound states at k = in as the dimension of 
Ker[J(m)T]. Since that is also equal to the dimension of Ker[J(i«)], we can say that the 
multiplicity of the bound state at k = in is given by the dimension of Ker[J(m)]. Let us 
use m K to denote the multiplicity of the bound state at k = in. We thus have 

m K = dimKer[J(m)]. (8.10) 
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Note that 1 < m K < n because the dimension of Ker[J(m)] cannot exceed n for the 
corresponding n x n matrix J(k). 

We summarize the above observations on bound states in the following theorem. 

Theorem 8.1 Consider the self adjoint matrix Schrddinger operator with the self adjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Let 
f(k,x), <p(k,x), and J(k) be the corresponding Jost solution, the regular solution, and the 
Jost matrix, appearing in (2.1), (2.5), and (3.1), respectively. Then: 

(a) We have a bound state at k = in for some positive k if and only if Ker[J(m)] is 
nontrivial or equivalently if and only if det[J(iK)] = 0. 

(b) The multiplicity m K of the bound state at k = in is finite, and in fact it is equal to 
the dimension of Ker[J(in)]. 

(c) A bound-state column-vector solution to (1.1) at k = in must be equal to f(in, x) (5 for 
some nonzero column vector (3 G Ker[J(in)^]. Similarly, a bound-state column-vector 
solution to (1.1) at k = in must be equal to <p(iK, x) a for some nonzero column vector 
a G Ker[J(iK,)]. 

(d) If k = ik corresponds to a bound state, then there is a bijection a ^ 13 between 
Ker[J(iK,)] and Ker[J(iK,y] in such a way that (p(in,x) a = f(iK, x) /3. 

We will next analyze the behaviors of the Jost matrix J(k) and of its inverse at a 
bound state k = in. One of our goals is to prove that the multiplicity m K of the bound 
state is equal to the multiplicity of the zero of det J(k) at k = in. We will use an overdot 
to indicate the derivative with respect to k. 

Proposition 8.2 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Let 
f(k,x) and J(k) be the corresponding Jost solution and the Jost matrix, appearing in (2.1) 
and (3.2), respectively. Assume that there is a bound state at k = in for some positive k. 
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Then, for each fixed x G R + we have 



f(lK, x)\ 



df(-k*, x y 




(8.11) 



dk 



df'(-k*,x)i 




/"(•K,x)t, 



(8.12) 



dk 



= f'{in) ] A - f(in, 0) + S. 



(8.13) 



PROOF: As stated in Section 2, f(k, x) is analytic in k G C + for each fixed x G R + . Thus, 
we have the Taylor series expansion 



Replacing k by — k* in (8.14) and by taking the adjoint of both sides of the resulting 
expansion, we get 



f(-k*,x)^ = f(iK,x)i - (k-iK)f(iK,x)i +0((k-i K ) 2 )), k^iK. (8.15) 



The first and second terms in the expansion on the right-hand side in (8.15) yield the 
equalities in (8.11). The equalities in (8.12) are established in a similar manner by exploit- 
ing the analyticity of f'(k,x) and f"(k,x) in k G C + for each fixed x G R + . By taking 
the /c-derivative of both sides of (3.2) and using the second equality in (8.11) and the first 
equality in (8.12), we obtain (8.13). | 

Theorem 8.3 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6) -(1.8) and with the potential V satisfying (1.3) and (1.5). Let 
f(k, x), <p(k,x), and J(k) be the corresponding Jost solution, the regular solution, and the 
Jost matrix, appearing in (2.1), (2.5), and (3.1), respectively. Assume that there is a bound 
state at k = in for some positive k. For any constant column vector a G Ker[J(iK)], let f3 
be the corresponding unique constant column vector indicated in Theorem 8.1(d). Then: 



f(k 7 x) = f(in, x) + (k — in) f(in, x) + 0((k — in) 2 )), k — > in. 



(8.14) 




(8.16) 
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and hence ft J(in)a ^ unless a = 0. 

PROOF: The Jost solution f(k,x) satisfies (1.1) and hence 

f"(k, x) + k 2 f(K x) = V(x) f(k, x). (8.17) 

By taking the /c-derivative, from (8.17) we get 

f"(k, x) + k 2 f(k, x) + 2kf(k, x) = V(x) f(k, x), (8.18) 

and by replacing k by — k* in (8.17) and then taking the adjoint we get 

f"(-k* , x)t + k 2 f(-k\ x)t = f(-k*,x)W(x), (8.19) 

where we have used the selfadjointness of V given in (1.3). Evaluating (8.18) and (8.19) 
at k = in, with the help of (8.11), we get 

f"(lK, X) - K 2 f(iK, X) + 2lKf(lK, x) = V(x) f(m, x), (8.20) 

f"(iK, X) j - K 2 f(iK, X) j = f(iK, x) j V(x). (8.21) 

Premultiplying (8.20) by f(iK 7 xY and postmultiplying (8.21) by f{iK,x) and taking the 
difference of the resulting equations, we obtain 

= -2w/(w, x)V(i«, x). (8.22) 

Premultiplying (8.22) by ft and postmultiplying it by (3, we integrate the resulting equation 
over x E R + . We then get 

poo 

-ftf(m, 0)V'(iK, 0)0 + ftf'{iK, 0) t /(z«, 0)0 = -2ik / dx [f(in, x)^[f(in, x)(3], 

Jo 

(8.23) 

where we have used Theorem 8.1(c) with the fact that f(w, x)j3 is a bound-state column- 
vector solution and hence it is square integrable, and we have also used the fact that 
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a 

d,x 



f(iK, x) ] f'(in, x) - f'(iK, x) ] f(in, x) 



the quantity inside the brackets in (8.22) vanishes as x — > +00. The latter property is a 
consequence of the exponential decay to zero of f(iK,x)/3 and of f'(iK,x)/3 and can be 
established with the help of (2.1) at k = in. Multiplying (8.23) on both sides by i and 
using (8.6) we get 

— i a^cp(iK, 0) ^f'{iK, 0)f3 — a^cp'(iK, 0) ^f(iK, 0)/?J = 2k J dx[(p(iK,x)a]^[(p(iK,x)a]. 

(8.24) 

Using (2.5) and by taking the adjoint of both sides in (8.24) we get 

r°° 

ift f'(iK,0)^A- fiiK^B a = 2n dx[(p(iK,x)a]^[(p(iK,x)a]. (8.25) 
L J Jo 

Comparing the left-hand side of (8.25) with (8.13), we obtain (8.16). Finally, since k > 

we see that the right-hand side in (8.16) is positive if a 7^ and is equal to zero if a = 0. 

Thus, from the left-hand side of (8.16) we conclude that ^ J(in)a = only when a = 0. | 

Theorem 8.4 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6) -(1.8) and with the potential V satisfying (1.3) and (1.5). Let 
f(k,x), <p(k,x), and J(k) be the corresponding Jost solution, the regular solution, and the 
Jost matrix, appearing in (2.1), (2.5), and (3.1), respectively. Assume that there is a bound 
state at k = in for some positive k. Then, J(k)~ 1 has a simple pole at k = in. 

PROOF: By Theorem 8.1(a), the determinant of J(k) vanishes at k = m, and hence 
Theorem 3.1(a) implies that J(k)~ l is analytic in a deleted neighborhood of k = in with a 
pole of some finite order p at k = in. If the pole at k = in were not simple, then for p > 2, 
in some neighborhood of k = ik we would have the expansions 

J(k) = J(ik) + (k- in) + 0({k - m) 2 ), (8.26) 

JW - 1 = (^ + (^i^ + ''' + ^ + ^ + ^ (8 ' 27 » 

Using (8.26) and (8.27) in J(/c)J(/c) _1 = 7 n , we would obtain 

J(ik) N- p = 0, J(i«) N- p+1 + N- p = 0. (8.28) 
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From the first equation in (8.28) we see that each column of N_ p would have to belong to 
Ker[J(m)]. For each nonzero column of iV_ p , by denoting that nonzero column with a as 
in Theorem 8.3, from (8.28) we would get the column-vector equation 

J(ik)( + J(iK)a = Q, (8.29) 

where a G Ker[J(m)] and ( is some column vector in C n .. Let (3 G Ker[J(w)T] be the 
unique column vector corresponding to a as stated in Theorem 8.3. Thus we would have 
J(inY (3 = or equivalently 

ft J(in) = 0. (8.30) 
Let us premultiply (8.29) by ft and use (8.30) in order to obtain 

ftj(in)a = 0. (8.31) 

Using Theorem 8.3 in (8.31) we see that we must have a = and hence N- p = for p > 2. 
Thus, from (8.27) we conclude that J{k)~ x must have a simple pole at k = in. | 

Having established that the expansion (8.27) contains only a simple pole as 

J(fc)" 1 = N ~ 1 + N + (k- in) Ni + 0((k - in) 2 ), k m, (8.32) 
k — in 

we would like investigate the term N-\ further. One of our goals is to relate the multiplicity 
of the bound state at k = in to the multiplicity of the zero of det J(k) at k = in and to 
show that those two multiplicities are equal to each other. Recall that the multiplicity 
m K of the bound state at k = in is defined as the number of linearly independent column 
vectors that are square-integrable column-vector solutions to (1.1) at k = in and that 
also satisfy the boundary condition (1.6)-(1.8). From (8.10) we see that m K is equal to 
the dimension of the kernel of J (in). Our goal is to prove that m K is also equal to the 
multiplicity of the zero of det J(k) at k = in. 

Theorem 8.5 Consider the selfadjoint matrix Schrddinger operator with the selfadjoint 
boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). Let 

37 



J(k) be the corresponding Jost matrix appearing in (3.1). Assume that there is a bound 
state at k = in for some positive k. Then, we have 



detJ(fc) = c 3 (k -m) mK [l + Q(k -in)], 



k — > in, 



(8.33) 



where C3 is a nonzero constant and m K is the positive integer appearing in (8.10) and 
denoting the multiplicity of the bound state at k = in. Consequently, the order of the zero 
of detJ(k) at k = in is equal to m K . 

PROOF: From (8.10) we know that the geometric multiplicity of the zero-eigenvalue of 
J(ik) is equal to m K . We will proceed as in Section 6 of [2], and hence we will omit some 
of the details by referring the reader to [2]. Using a similarity transformation J (in) h-> 
S^ 1 J(ik)Si with an appropriate invertible matrix <Si, we will transform J (in) to a Jordan 
canonical form. Let us assume that there are v K Jordan chains and hence the Jordan 
canonical form of J (in) contains v K Jordan blocks. Let us use A s to denote the eigenvalue 
of J (in) associated with the sth Jordan chain, where we realize that the eigenvalues may 
be repeated and hence there may be more than one Jordan block for a given eigenvalue A s . 
Let us use J ns (X s ) to denote the sth Jordan block, where we assume that the matrix size 
of that block is n s xn s . Without loss of generality we can assume that the first m K Jordan 
chains all belong to the zero eigenvalue of J {in). Let us use fi K to denote the algebraic 
multiplicity of the zero eigenvalue of J(m). Thus, we assume that the number of nonzero 
eigenvalues (including multiplicities) of J (in) is n — \i K . As a result, the first m K Jordan 
blocks each have the form 



ro 1 



01 



1 







8=1,..., 



rn 



(8.34) 







1 











and the remaining Jordan blocks associated with the nonzero eigenvalues of J {in) have 
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the form 



X s 1 .. 

A s 1 .. 

0.. 

L 0.. 








A s 1 



s = m K + 1, 



X a . 

with nonzero diagonal entries A s . The Jordan canonical form of J (in) is then given by 

5r x J(i«) 5i = ©^ 1 J ns (A s ). 

Next, let us move all the entries with 1 appearing in the superdiagonal in the first m K 
Jordan blocks in (8.34) and collect those entries into the — m K ) x {p> K — m K ) identity 
matrix I flK _ rtlK . This can be achieved by using the matrices P4 and P5 given by 





"n 4 




"n 5 






, P 5 = 






. In-» K _ 


. /n- MK . 



for some permutation matrices II4 and II5 that affect only the first fi K columns and fi K rows, 
respectively, of the matrices on which they operate. The combined matrix transformation 
J (in) 1 — y P^S^ 1 J {ik)S\P^ results in the upper-triangular matrix given by 

P 5 S^ J(m)<SiP 4 = diag{0 mK , Ip K -m K , Jn^ K+1 (A Mk+ i), . . . , J n „ K (K K )}, (8.35) 

where we recall that mK denotes the m K x m K zero matrix. Let us define the (n — /i K ) x 
(n — fx K ) matrix do as 



d := diag{/ MK _ mK , J nfiK+1 (A Mlt+ i), . . . , J n „ K (A V|t )}. 



5.36) 



The matrix do is invertible because it is an upper-triangular matrix with nonzero diagonal 
entries. Using (8.36) in (8.35) we obtain the block decomposition 



P 5 S 1 J(ik)SiP4 = diag{0 mK ,d o }. 



(8.37) 
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Comparing (8.26) and (8.37) we see that 
P 5 S^ 1 J(k)S 1 P 4 = diag{O mK , d } + (k-iK) 
where we have let 



ai bi 
ci di 



+ 0((/c-m) 2 ), fc-H/c, (8.38) 



:= P 5 S^ j(iK)SiP 4 . 



ai bi 
ci di 

From Theorem 8.4 we know that J(/c) _1 has a simple pole at k = in and hence with the 
help of (8.32) we get 



(Pssr 1 ms.p,)- 1 = j-L- 



where we have defined 

ni n 2 
n 3 n 4 



ni n 2 




mi 


m 2 




+ 






n 3 n 4 




m 3 


m 1 




m. 


m 2 


-1S1-F5 ' 










m 1 



+ 0(k-iK), k ->■ i«, (8.39) 



with iV_i and iVo being the matrices appearing in (8.32), with some m K xm K block matrices 
ni and mi, some (n — m K ) x (n — m K ) block matrices 114 and 1114, and the remaining block 
matrices of appropriate sizes. Using (8.38) and (8.39) in the matrix identities 

(P 5 5r 1 J(A;)5iP4)" 1 {P 5 S^ 1 J(k)S 1 P 4 ) = /„, 

(P 5 5r 1 J(A;)5iP4) {P 5 S^ 1 J(k)S 1 P 4 )~ 1 =I n , 
we obtain 

n 1 no n 1 no 

(8.40) 



(8.41) 
(8.42) 



ni 


n 2 


diag{0 mK 


,d } 


= 


n , diag{0 mK ,d } 


ni 


n 2 




n 3 


n 4 
















n 3 


n 4 






ni 


n 2 




ai 


bi" 


+ 


mi m 2 


diagjO^ 




Ad 




n 3 


n 4 




_ c i 


di. 




m :j m 1 














ai 


bi" 




ni 


n 2 


+ diag{0 mre ,d } 


mi 


m 2 








Cl 


di 




n 3 


n 4 








m 3 


m 4 







Because d is invertible, from (8.40) we see that 



n 2 = 0, n 3 = 0, n 4 = 0, 
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(8.43) 



for some zero matrices of appropriate sizes. Using (8.43) in (8.41) and (8.42) we get 

niai = I mK , nibi + m 2 d = 0, m 4 d = I n -m K , Cini + d m 3 = 0, 
which establishes the invertibility of the block matrix ai and also implies 

n^af 1 , m 2 = -af^ido 1 , m 4 = dg \ m 3 = -dg 1 CiaJ" 1 . (8.44) 



Using (8.44) in (8.39) we obtain the expansion 
(PsSf 1 J(/c)<SiP4) _1 = 



aiVm K +0(k-iK)] , ! 

— a 1 diUq 



— m 

-do 1 c i a^ 1 d^ 1 



(8.45) 



From (8.45) we see that 

det ((PscSr 1 J(A;)5 1 P 4 )" 1 ) = J ^ - t 1 + 0(k - ik)} , k ->■ zk, 
or equivalently 

det J(fc) = ^^pt ! (fc - *«) m " [1 + ~ »«)] » fc->i/c. (8.46) 

aet (P4P5) 

Thus, (8.46) establishes (8.33) with c 3 given by 

c 3 := (detP 4 )(detP 5 )(detai)(detd ), 

which is nonzero due to the fact that P 4 and P5 are some n x n permutation matrices and 
hence their determinants are either 1 or —1, and the matrices ai and do are invertible and 
hence their determinants are nonzero. | 

We remark the similarity between Theorem 8.5 and Corollary 6.2 and the similarity 
between (8.33) and (6.5). 

Let us note that the transformation specified in Proposition 4.1(a) on the boundary 
parameters A and B does not affect the boundary condition (1.6)-(1.8). This is because 
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(A, B) h-> (AT, BT) for an invertible matrix T results in a premultiplication of both sides 
of (1.6) by Tt as well as a premultiplication by and a postmultiplication by T of both 
sides of (1.7) and (1.8). Thus, as seen from (3.2), the potential V and the boundary 
parameters A and B cannot uniquely determine the Jost matrix J(k), but they determine 
J(k) uniquely up to a postmultiplication by an invertible matrix T. However, such a 
nonuniqueness does not affect the zeros in C + of the determinant of J(k) because det J(k) 
and det [J(k)T] have the same set of zeros. Hence, the bound states are not affected by 
such a nonuniqueness, and the bound states are uniquely determined by the potential V 
and the boundary parameters A and B appearing in (1.6)-(1.8). 

The following result is relevant in establishing the finiteness of the number of bound 
states. 

Theorem 8.6 Consider the Schrddinger operator corresponding to (1.1) with the self ad- 
joint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). 
Let J(k) be the corresponding Jost matrix defined in (3.1). Then, the only zeros of det J (k) 
in C+ \ {0} can occur on the positive imaginary axis, and the number of such zeros, which 
we denote by N (without counting multiplicities), is finite. 

PROOF: By Theorem 3.1(a) we know that J(k) is analytic in C + and continuous in 
C+. Thus, det J(k) possesses the same properties. Because of the selfadjointness of the 
Schrddinger operator, the bound-state /c-values, i.e. the zeros of det J(k) in C+ \ {0}, can 
occur either on the real axis or on the positive imaginary axis. By Proposition 3.1(c), J(k) 
is invertible for k G R\ {0} and hence those zeros can only occur on the positive imaginary 
axis. Let us use H to denote the set of zeros of det J(k) on the positive imaginary axis. 
Because of (7.17), H is a bounded set. Furthermore, (6.5) implies that det J(i() ^ for 
< C < kq for some positive Ko-value. Thus, it follows that H C [ikq, ib] for some positive 
b. We must prove that H is a finite set. If it were not a finite set, being bounded, H would 
have to have an accumulation point in [ikq, ib]. However, the analyticity of det J(k) in C + 
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would then require det J(k) = in C + , contradicting (7.17). | 

Let us assume that the A" distinct zeros of detJ(k) on the positive imaginary axis 
occur at k = iKj with j = 1, . . . , N. If there are no bound states, then we have N = 0. If 
there any bound states, as stated in Theorem 8.6, the positive integer A" is finite. Let us 
use m Kj to denote the multiplicity of the bound state at k = iKj. As in (8.10), we have 

m Kj = dimKer[J(mj)], 

and hence m Kj is a positive integer not exceeding n. From Theorems 8.1 and 8.6 we conclude 
that, the number of bound states including the multiplicities, A/", is a finite number and 
given by 

N 

i=i 

From Theorem 8.5 it follows that the multiplicity of the zero of det J(k) at k = iKj is the 
same as the multiplicity m Kj of the bound state at k = iKj. Thus, from (8.47) we have the 
following result. 

Corollary 8.7 Consider the Schrddinger operator corresponding to (1.1) with the self ad- 
joint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). 
Let J(k) be the corresponding J ost function given in (3.1), and let M be the corresponding 
number of bound states (including multiplicities), as indicated in (8.47). Then, Af is also 
equal to the number of zeros (including multiplicities) of det J (k) in C + . 

In Section 9 we will relate Af to the change in the argument of the determinant of the 
scattering matrix S(k) along the positive real axis. 

9. LEVINSON'S THEOREM 

In this section we establish Levinson's theorem for the selfadjoint matrix Schrddinger 
operator with the selfadjoint boundary condition (1.6)-(1.8) and with the potential V 
satisfying (1.3) and (1.5). We do this by relating the argument of the determinant of the 
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scattering matrix S(k) defined in (3.8) to the number (including multiplicities) of bound 
states M given in (8.47). We achieve our goal by applying the argument principle to the 
determinant of the Jost function J(k) given in (3.1). 

Let us define h(k) as 



The region we will use in the argument principle is the region whose boundary is C 6j r, 
which consists of four pieces as given by 



The first piece (-R, — e) is the directed line segment on the real axis for some small positive 
e and for a large positive R, with the direction of the path from —R + iO to — e + iO. The 
second piece C e consists of the upper semicircle centered at the origin with radius e and 
traversed from the point — e + iO to the point e + zO. The third piece (e, R) is the directed 
line segment of the positive real axis from e + iO to R + iO. The fourth piece Cr is the upper 
semicircle centered at the origin with radius R and traversed from the point R + iO to the 
point — R + iO. The analyticity of h(k) in our region and its continuity in the closure of our 
region follows from Theorem 3.1(a). By choosing R large enough and by choosing e small 
enough, from Theorem 8.6 we know that the only zeros of h(k) in our region can occur on 
the positive imaginary axis at N distinct points k = iKj for some nonnegative integer N 
and that h(k) does not vanish on the boundary of our region. 

Let us use arg[7i(fc)]| c for the change in the argument of h(k) along a path C, and let 
us recall that an overdot indicates the /c-derivative. 

Proposition 9.1 Consider the Schrddinger operator corresponding to (1.1) with the self- 
adjoint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). 
Let J(k) and S(k) be the corresponding Jost and scattering matrices defined in (3.1) and 
(3.8), respectively. Then, the change in the argument of det S(k) along the directed path 



h(k) := det J(k). 



(9.1) 



C e , R := {-R, -e) UC £ U (e, R) U C R . 



(9.2) 
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(e, R) and the change in the argument of det J(k) along the directed paths (-R, — e) and 
(e, R) are related to each other as 

arg[det S (k)} \ {e R) = -arg[det J (k)} \ (e R) - arg[det J (k)] | _ e) • (9.3) 

PROOF: From (3.8) we have 

det^fc) = (-l)- d ^ t J j~^ , keR\{0}. (9.4) 
From (3.11) we get | det5'(/c)| = 1 and hence 

detS(fc) = e iais[dets{k)] , k G (e,R), (9.5) 

| det J(-k)\ = | det J(k)\, keH\{0}, (9.6) 

where we have also used (9.4) to obtain (9.6). By Theorem 3.1(c) we have | det J(k)\ ^ 
for k E R \ {0}, and hence with the help of (9.6) we get 

det J{k) = | det J(fc)|e iarg[detJ(fc)1 , ke(e,R), (9.7) 

det J(-fc) = | det J(A;)|e iarg[detJ( - fc)1 , ke{e,R). (9.8) 
Using (9.5), (9.7), and (9.8) in (9.4) we get 

g iarg[det S(k)] _ ^_^n e i arg[det J(-fc)] g -i arg[det J(k)] k E {t R) 

from which we obtain (9.3). | 

In the following proposition, we provide the change in the argument of h(k) along the 
pieces of paths appearing in (9.2). 

Proposition 9.2 Consider the Schrddinger operator corresponding to (1.1) with the self- 
adjoint boundary condition (1.6)-(1.8) and with the potential V satisfying (1.3) and (1.5). 
Then, the function h(k) defined in (9.1) satisfies 

lim lim I dk = 2niM, (9.9) 
e _> + R^+oo J Ce R h(k) 
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f h(k) 

lim / dk — — = iri(n M + n N ), (9.10) 
r ^ +00 J Cr h(k) 

lim / dk = -irifi, (9.11) 

V(-fl,-e)U(e,fl) A H'V V y " 

where Af is the nonnegative integer appearing in (8.47), the paths C e ,R, C e , and Cr are 
those in (9.2), %/ and tin are the nonnegative integers defined after (4-2), and fx is the 
(algebraic and geometric) multiplicity of the eigenvalue +1 of the zero-energy scattering 
matrix S(0), with S(k) being the scattering matrix defined in (3.8). 

PROOF: Because of Corollary 8.7, the number of bound states Af (including multiplicities) 
is equal to the number of zeros of det J(k) (including multiplicities) in C + . We get (9.9) by- 
applying the argument principle to h(k) along the closed path C €j r and by using the fact 
that Af is the number of zeros (including multiplicities) of h(k) inside C ej R. We note that 
(9.10) directly follows from (7.17), and (9.11) directly follows from (6.5). Finally, (9.12) is 
obtained with the help of (9.6). | 

We next state Levinson's theorem. 

Theorem 9.3 Consider the Schrddinger operator corresponding to (1.1) with the self ad- 
joint boundary condition (1.6) -(1.8) and with the potential V satisfying (1.3) and (1.5). 
The number Af of bound states (including multiplicities) appearing in (8.47) is related to 
the argument of the determinant of the scattering matrix S(k) defined in (3.8) as 

arg[det S(0 + )] - arg[det 5(+oo)] = n (2Af + n - n M - n N ) , (9.13) 

where fx is the (algebraic and geometric) multiplicity of the eigenvalue +1 of the zero-energy 
scattering matrix S(0), and um and are the nonnegative integers defined after (4-2). 

PROOF: By Proposition 3.3, the determinant of S{k) is continuous on R and hence the 
left-hand side in (9.3) is given by 

lim lim (argfdet S(k)\ L m ) = -argfdet S(0 + )] + argfdet Sf+oo)], (9.14) 
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By combining the results in (9.10)-(9.12), and by using the relationships among (9.3), 
(9.12), and (9.14), we evaluate the sum of the integrals in (9.10)-(9.12) in the limit as 
e — > + and R — > +oo. That sum then must be equal to the value of the integral given in 
(9.9), resulting in (9.13). | 

Let us comment on (9.13) and how it is related to Levinson's theorem appearing in 
the literature elsewhere. In our analysis of the selfadjoint Schrodinger operator with the 
general selfadjoint boundary condition (1.6)-(1.8), the unperturbed Hamiltonian is chosen 
to satisfy the Neumann boundary condition. Such a choice is compatible with the time- 
dependent derivation of the scattering matrix and is motivated [14] by applications in 
quantum wires. One consequence of that choice is apparent in the large-/c limits of S(k). 
As seen from (5.7), (5.10) and (7.20), we have S(k) = S(oo) + 0(l/k) as k ->• ±oo, with 

S(oo) = M diag{/ nM , -I nD , I nN } M+ , (9.15) 

where M is the unitary matrix appearing in (4.11). As a result of (9.15), the argument of 
the determinant of S(k) as k — > +oo is given by 

lim arg[detS(fc)] = argfdet S(oo)] = (-l) n °7r + 27rj, j = 0, ±1, ±2, . . . . 

k— >-+oo 

One can choose a branch with arg[det S'(oo)] = if no is even and a branch with 
arg[det S'(oo)] = tt if no is odd. So, in the purely Dirichlet case, i.e. when tim = ^jv = 
and n_o = n, from (9.15) we get S(k) = —I n + 0(1/ k) as k — > ±oo. On the other 
hand, in the literature [1,21,22] dealing solely with the Dirichlet case, it is customary 
to choose the unperturbed Hamiltonian to satisfy the Dirichlet boundary condition, yield- 
ing S(k) = I n + 0(1/ k) as k — > ±oo. Hence, in the literature dealing solely with the 
Dirichlet case, it is also customary to use the particular branch of the argument function 
for det S(k) in such a way that that argument takes the zero value at k = +oo. In fact, in 
such a case, it is customary to let 

detS(fc) = e 2i5s{k \ fce(0,+oo), 
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with 5s (+00) = 0. Then, Levinson's theorem in such a case is given by 

5 5 (0+) = 7r(V+|). (9.16) 

In particular, in the scalar case , we have (9.16) with \i = 1 in the exceptional case and 
fj, = in the generic case [21]. 

Acknowledgment. The research leading to this article was supported in part by Consejo 
Nacional de Ciencia y Tecnologia (CONACYT) under project CB2008-99100-F and by the 
Department of Defense under grant number DOD-BC063989. 

References 

[1] Z. S. Agranovich and V. A. Marchenko, The inverse problem of scattering theory, 
Gordon and Breach, New York, 1963. 

[2] T. Aktosun, M. Klaus, and R. Weder, Small-energy analysis for the self-adjoint matrix 
Schrddinger equation on the half line, J. Math. Phys. 52, 102101 (2011). 

[3] T. Aktosun and R. Weder, Inverse spectral- scattering problem with two sets of discrete 
spectra for the radial Schrddinger equation, Inverse Problems 22, 89-114 (2006). 

[4] G. Berkolaiko, R. Carlson, S. A. Fulling, and P. Kuchment (eds.), Quantum graphs and 
their applications, Contemporary Mathematics, 415, Amer. Math. Soc, Providence, 
RI, 2006. 

[5] J. Boman and P. Kurasov, Symmetries of quantum graphs and the inverse scattering 
problem, Adv. Appl. Math. 35, 58-70 (2005). 

[6] P. Exner, J. P. Keating, P. Kuchment, T. Sunada, and A. Teplyaev (eds.), Analysis on 
graphs and its applications, Proc. Symposia in Pure Mathematics, 77, Amer. Math. 
Soc, Providence, RI, 2008. 

[7] N. I. Gerasimenko, The inverse scattering problem on a noncompact graph, Theoret. 

48 



Math. Phys. 75, 460-470 (1988). 

[8] N. I. Gerasimenko and B. S. Pavlov, A scattering problem on noncompact graphs, 
Theoret. Math. Phys. 74, 230-240 (1988). 

[9] B. Gutkin and U. Smilansky, Can one hear the shape of a graph? J. Phys. A 34, 
6061-6068 (2001). 

[10] M. S. Harmer, Inverse scattering for the matrix Schrddinger operator and Schrddinger 
operator on graphs with general self-adjoint boundary conditions, ANZIAM J. 44, 161- 
168 (2002). 

[11] M. S. Harmer, The matrix Schrddinger operator and Schrddinger operator on graphs, 
Ph.D. thesis, University of Auckland, New Zealand, 2004. 

[12] M. Harmer, Inverse scattering on matrices with boundary conditions, J. Phys. A 38, 
4875-4885 (2005). 

[13] V. Kostrykin and R. Schrader, Kirchhoff's rule for quantum wires, J. Phys. A 32, 
595-630 (1999). 

[14] V. Kostrykin and R. Schrader, Kirchhoff's rule for quantum wires. II: The inverse 
problem with possible applications to quantum computers, Fortschr. Phys. 48, 703-716 
(2000). 

[15] P. Kuchment, Quantum graphs. I. Some basic structures, Waves Random Media 14, 
S107-S128 (2004). 

[16] P. Kuchment, Quantum graphs. II. Some spectral properties of quantum and combi- 
natorial graphs, J. Phys. A 38, 4887-4900 (2005). 

[17] P. Kurasov and M. Nowaczyk, Inverse spectral problem for quantum graphs, J. Phys. 
A 38, 4901-4915 (2005). 

[18] P. Kurasov and F. Stenberg, On the inverse scattering problem on branching graphs, 

49 



J. Phys. A 35, 101-121 (2002). 

[19] B. M. Levitan, Inverse Sturm- Liouville problems, VNU Science Press, Utrecht, 1987. 

[20] V. A. Marchenko, Sturm- Liouville operators and applications, Birkhauser, Basel, 1986. 

[21] R. G. Newton, Scattering theory of waves and particles, 2nd ed., Springer, New York, 
1982. 

[22] R. G. Newton and R. Jost, The construction of potentials from the S -matrix for 
systems of differential equations, Nuovo Cimento 1, 590-622 (1955). 



50 



